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مه ‌ 

هنگامی که ابعاد مواد تا ثانومت رکاهش می‌پابد بر حصوصیات ذره مثل نقطه ذوب و در نتیجه نمودار فازی تعادلی اث رگذار است. در این مطالعه برنامه‌ای 
در محیط نرم‌افزاری متلب تدوین شده است. با اعمال روش‌های مبتنی بر ترمودینامیک سطح,. از این برنامه برای رسم نمودارهای فازی نانوذرات, 
برر سی تاثیر شعاع نانوذرات بر نمودار فازی تعادل یآنها و مقایسه آن با حالت توده استفاده شده است. با توجه به نتایج حاصل استنباط می شود که با 
کاهش اندازه ذرات: حطوط جامد و مذاب در نمودار فازی تعادلی به طرف پایین حرکت م یکنند. علاوه بر این حلالیت حالت جامد در مقیلس ثانو تا 


حد زیادی افزایش يافته و مناطق دو فازی باریک‌تر می‌شوند. 


واژه های کلیدی نانوذرات اندازه ذره نقطه ذوب. نمودار فازی. نرم‌افزار متلب. 
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مقد مه 

بکارگیری مواد نانو در سال‌های اخیر توسعه چشم‌گیری 
یافته‌است که از آن جمله می‌توان به صنایع مختلفی از جمله 
الکترونیک. کاتالیز سرامیک. ذخیره‌سازی مغناطیسی 
اطلاعات. اجزای ساختمانی و ... اشاره نمود؛ به عنوان مثال 
مینیاتوری نمودن تجهیزات الکترونیکی نیاز به جایگزین 
کردن قطعات با ابعاد در حد نانومتر دارد. با کاهش ابعاد تا 
این حدود. مواد خواص مکانیکی و فیزیکی جالب و غیر 
منتظره‌ای از خود نشان می‌دهند که افزایش استحکام 
مکانیکی. افزايش ضریب نفوذ [1-3 بالا رفتن گرمای ویژه 
و مقاومت الکتریکی [4] ,[1] و کاهش نقطه ذوب [5] و در 
نتیجه تغییر نمودار فازی تعادلی از آن جمله‌اند. 

به‌منظور بررسی اثر اندازه بر نقطه ذوب نانوذرات و 
نمودار فازی تعادلی آنها مطالعات بسیاری صورت گرفته 
است. به‌عنوان مثال در سال ۲۰۰۸ مدلی برای پیش‌بینی نقطه 
ذوب نانوذرات بر اساس عدد کوئوردیناسیون متوسط ارائه 
شده است که نشان می دهد با توجه به اينکه دمای ذوب 
متناسب با عدد کوئوردینانسیون می‌باشد و این عدد در 
نانوذارت نسبت به حالت توده کم‌تر است. در نتیجه دمای 
ذوب نانوذرات با کاهش شعاع ذره روند نزولی را نشان می- 
دهد [5]. لی و همکارانش اثر اندازه ذرات بر روی دمای ذوب 
اس ۸۵0 را پررسی کردند. نتایج مطالعات آنها نشان 
می‌دهد که قطر ذرات با دمای ذوب ارتباط مستقیم و خطی 
دارد [6]. 

در تحقیق دیگری که توسط گارزل و همکارانش 
بمنظور بررسی تاثیر اندازه و شکل ذرات بر روی نمودار 
فازی دوتایی نانوذرات ۸-۷ صورت گرفته. گزارش شده 
ایشا که انلازه درا هاش کیک اف سا کرد 
و همچنین ترکیب یوتکتیک نیز به مقادیر کمتر ) در نمودار 
فازی انتقال پیدا می کند [7]. مطالعات تاناکا بر روی اثر اندازه 
ذرات بر سیستم های آلیاژی دوتایی نشان می‌دهد که در اثر 
کاهش قطر ذرات علاوه بر کاهش دمای ذوب. منطقه حلالیت 
در نمودارهای فازی افزايش می یابد [8]. نتایج تحقیقات 
یانگ و همکارانش در پی چگونگی محاسبه نمودارهای فازی 


رسم نمودارهای فازی ثانو برای سیستم‌های دوتایی در محیط نرم افزاری متلب 


نانوذرات نیز منجر به یافتن رابطه خطی بین آنتالپی ذوب؛ 
دمای ذوب و انرژی بر هم کنش اتمی (462 با انرژی پیوستگی 
(ع۳) شد. آنها ادعا کردند با کاهش قطر ذرات. خحطوط جامد 
و خط مذاب در دماهای پایین‌تری نسبت به حالت بالک ماده 
قرار می‌گیرد و منطقه دو فازی نیز کاهش می‌یابد [9]. لیانگ 
و هایمینگ در تحقیقات خود به نتایجی مشابه نتایج یانگ 
دست یافتند. به اين ترتیب در این مدل‌ها می‌توان با 
جایگذاری آنتالپی ذوب. دمای ذوب و انرژی بر هم کنش 
اتمی نانوذره به جای مقادیر مربوط به توده آن نمودار فازی 
تعادلی نانوذرات را به دست آورد [11] ,[10]. 

هماتطوی هنن ی ای شا نا تاکن متل فا 
ش ی سرام نی تیوفان قاری تانید رات ار فیه اش 
هدف از این پژوهش ارائه الگوریتمی نوین به منظور رسم 
نمودارهای فازی نانو در سیستم های باقاعده و بی قاعده 
است که با پیاده‌سازی آن در نرم‌افزار متلب بتوان نمودارهای 
فازی را در حالت توده و نانو با دقت قابل قبولی در حداقل 
زمان ممکن رسم کرده و تاثیر معادلات ارائه شده در این 


مدل‌ها را بر نمودار فازی نانوذرات مشاهده نمود. 


موی ووشن ین 
همانطور که گفته شد هدف اصلی در این تحفیق علاوه بر 
پروشی یز تاه دانه بر تموفاوهاق قاری تمادی) جافت 
راه‌حلی است تا با حداقل هزینه و زمان بتوان نمودارهای 
تعادلی را در شرایطی به دور از آزمایشات پیچیده رسم کرد. 
در این راستا برنامه‌ای در محیط نرم‌افزاری متلب نوشته شده 
است. 
نرم افزار متلب (1۷۸1۸8) برنام‌ای کامپیوتری است 
که راغ کتاقی کانبات عزوت. ووهه: خه عظی 
سر و کار دارند» تهیه شده است. کار کردن با ماتریس‌ها در 
متلب بسیار ساده است. در حقیقت تمام داده‌ها در آن به شکل 
یک ماتریس ذخیره می‌شوند. 
در این مطالعه به منظور ساده سازی محاسبات. ابتدا 
معادلات انرژی آزاد گیبس تمام فازهای موجود بر حسب دما 


و درصد اتمی در نظر گرفته شد. سپس در هر دما نمودار هر 
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فازه بر حسب درصدهای جزئی ترسیم گردید. بررسی ها 
نشان داد نمودارهای فازهای مختلف. در برخی درصدهای 
ای شا کر نم وشات تخت انس ما مهو ننک 
از لحاظ ترمودینامیکی. فاز با انرژی آزاد گیبس کمتر پایداری 
بیشتری دارده در هر یک از این دماها. نمودار انرژی آزاد 
گیپس هر یک از فازهای موجود در تمامی درصدهای اتمی 
با یکدیگر مقایسه شد و فاز دارای حداقل انرژی به عنوان فاز 
پایدار در یک ماتریس تحت عنوان "ماتریس مینیمم" با ابعاد 
مشخص ذخیره گردید (جدول او شکل ۱). اين عمل برای 
تمامی دماها تکرار شد و در نهایت تمام منحنی‌ها به یک 
منحنی مینیمم تبدیل گردید. 

فز انامه خلت بل کروق زرط عمانن ی عی قظ از 


نقاط بعدی روی منحنی یک خط متصل می‌شود. سپس هر 
خط رسم شده به ازای تمامی مقادیر محور طولی (درصد 
اتمی) با تابع مینیمم مقایسه می‌شود.. تنها در صورتی که خط 
عبوری به طور کامل در تمام بازه‌ها پایین تر از تابع مینیمم 
قرار گیرد. خط مورد نظر یک مماس خواهد بود (شکل ۲). 
بدین ترتیب پس از پیدا کردن خط موردنظر. نقاطی از تابع 
مینیمم که خط مذکور برآن‌ها مماس گردیده است در یک 
ماتریس جمع‌آوری می‌گردد. در مرحله آخر نیز با رسم 
ماتریس نهایی (حاصل از نقاط مماس مشترک) نمودار فازی 
ال می توق 

در شکل (۳) فلوچارت و مراحل تر سیم نمودار فازی 


نمایش داده شده است. 


جدول ۱ پیدا کردن ماتریس مینیمم انرژی آزاد در یک دمای مشخص و تمام ترکیبات 
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شکل ۱ پیدا کردن منحنی انرژی آزاد مینیمم در یک دمای معین و در تمام بازه ترکیب شیمیایی 


سال سی و یکم» شماره دی ۱۳۹۹ 
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اس 


شکل ۲ محاسبه خطوط واصل در یک دمای معین و مقایسه آن با منحنی انرژی آزاد مینیمم در همان دما ( همانطور که در شکل پیداست تنها حط 1 


در تمام بازه پایین‌تر از تابع مینیمم قرار دارد؛ بنابراین نقاط 261و 262 از آن ذخیره می‌شود.) 


/ گرفتن مقادبر دماهای ذوب و / 


/ کام‌های مر بیط بد دما تر کیب / 


عه دبگر تقاط ماتریس 0 


1 


مقایسه معادلات خطوط واصل بدست آمذه نا مقادیر اثرژی وان 


همان دما 


ژ‌ ‌ 


محان یه ممادله خط واصل هر نقطله از ماتر یی میتمم شر هر دسا ۱ 


مینصم ذر هر دما و رطس آورمن خصو حط مماس بر متحنی آلرزرگ 


۲ 


٩‏ 1 رتیوت لت | محاسبه و ذخیره تفاط تلافی خطوط عماس با عنحتی انرژی آزاد 


ص .. میتیمم در هر دما 


شکل ۳ فلوچارت و مراحل رسم نمودار فازی در متلب 
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۵ دباگرامقازی وآقمی 
دیاگرام قازی محاسبه شده توسط برنامه 
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۵ دباگرام فازی وافعی 
دیاگرام فازی محاسبه شده توسط برنامه 
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شکل ۶ (الف) نمودار فازی سیستم 0-۷ (ب) نمودار فازی سیستم 1601-1061 


نتایج و بحث 
در این قسمت ابتدا تعدادی از نمودارهای فازی دوتایی توده 
به کمک نرم افزار رسم و بمنظور راستی آزمایی نرم‌افزار با 
نمودار فازی واقعی آنها مقایسه شده است. در ادامه نیز به 
کمک سه مدل ارائه شده در مقالات. نمودارهای فازی تعادلی 
نانو برای ذراتی با شعاع ۱۰ نانومتر ارائه و با حالت توده 


مقایسه شله. استخا: 


نتایج خروجی از نرم‌افزار در حالت توده 
در شکل (۶-لف و ب) نمودارهای فازی مربوط به دو 
سیستم (-00) و 1601-۷26 به کمک نرم‌افزار ارائه شده در 
این مقاله رسم شده است. همانطور که مشاهده می‌شود 
نمودارهای فازی که از حروجی نرم‌افزار بدست آمده تا حدود 
زیادی بر نمودار فازی واقعی این سیستم‌ها منطبق است. 
بررسی تاثیر اندازه ذره پر نمودار فازی 


سال سی و یکم» شماره دی ۱۳۹۹ 


در این بخش هدف. بررسی سه مدل پیش‌بینی نمودارهای 
فازی تعادلی نانوذرات. با اعمال فرضیات آنها در مقادیر 
ورودی نرم‌افزار می‌باشد. در مدل لیانگ فرضیاتی برای 
آتقالین ذوب. دمای ذوب و همچنین انرژی برهم‌کنش اتمی 
مطابق معادلات زیر در نظر گرفته شده است. 


1 0)وزن25 
)0 اس ت -] هه (م)آ-(ط) 1 


:,0(۳:۵( اه گأه(م)‎  ( 


٩2)(/۹ (۱-2 ۲/۳ ٩ - 010 (۳ 


در این روابط (0),رلگ ()یلط و ()یآ ‏ ()م و 
(200(446) به ترتیب آنتالبی ذوب. دمای ذوب و انرژی 
برهم‌کنش توده ماده و نانوذره مربوط به آن می‌باشد. علاوه 
بر این 1 ابت گازها. (0),ز,٩‏ آنتروپی ارتعاشی» 10 شعاع 


نانوذره و 100 نیز قطر بحرانی است که در آن تمامی اتم‌ها بر 
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سطح واقع می‌شوند و می‌توان از رابطه 2)3-0(18-(1 محاسبه 
می‌شود. به این ترتیب می‌توان با دادن مقادیر آنتالبی ذوب. 
دمای ذوب و انرژی برهم‌کنش اتمی مربوط به نانوذره به 
جای توده ماده به نرم افزان نمودار فازی نانوذارت را به دست 
آورد. مقادیر ذکرشده برای نانوذره آلیاژی ژرمانیم-سیلیسیم 
با قطر ۱۰ نانومتر و حالت توده آن در جدول (۲) آورده شده 
است و نمودار فازی در هر دوحالت به کمک نرم‌افزار در 
شکل (۵) رسم شده است. 

همانظزن که امه می‌شود با کاهقن آنذازه که ۶ 
محدوده نانو دمای ذوب کاهش می‌یابد. در نتیجه خحطوط 
جامد و مذاب به طرف پایین حرکت می‌کنند. هر چه میزان 
کاشین قظر مش قزر باقن سایانن ضظ این ی فیط بوارت 
نیز بیش‌تر اخواهد بود. علاوه بر این منطقه دو فازی کوچک- 
تر شده‌است. با توجه به رابطه (۳) با در نظر گرفتن دو برابر 
لایخ آنن بای قظر تا خی فقدای این مزر کرو 
محلول با قاعده به محلول ایده‌آل تبدیل می‌گردد. در اين 
صورت منطقه دوفازی دیگر در نمودار فازی مشاهده نمی- 
شود. از آنجایی که (),باط و (42/0 هر دو تابع اندازه ذره 
هستند» این مدل نسبت به سایر مدل‌ها کاهش بیش‌تری را در 
موقعیت خطوط جامد و مذاب و در اندازه منطقه دو فازی 


پیش‌بینی می‌کند [10]. 
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جدول ۲ اطلاعات مورد استفاده جهت ترسیم شکل (۳ [10] 
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در مدل هایمي نگ (عطنصعنه1؟) ارتباط بین آذتالپی 


ذوب. دمای ذوب و انرژی برهم‌کنش اتمی نانوذره و توده 


ماده به صورت زير تعریف شده است: 
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که در این رابطه ٩‏ آنترویی تبخیر می‌باشد [11]. 

مدل دیگری نیز توسط تاناکا ارائه شده است. در این 
مدل با فرض محلول باقاعده تاثیر ریز شدن ذرات بر نمودار 
حاوی محلول جامد بررسی شده است. با در نظر گرفتن فاز 
جامد خالص به عنوان حالت مرجع انرژی گیبس, برای یک 
سیستم آلیاژی با ذرات کوچک با شعاع 7 می‌توان انرژی 


گیبس کل را برای فازهای مذاب و جامد به صورت زیر 


نوشت: 
۱( ونلععه؟ ای ۳ هناگی ۳ ونالهاه کی 
2 ام5عع) ناگی ۳ امکعتاتگی ت امکهه کی 


در روابط فوق؛ تیان لع ۱,۱ و الستی انم انرژی های 


کشس وه فازهای مذاب و جامد بوده که در 00 < ۲ برابر 


با وناباهاه ی ۳ امکرلهاه ۸ خواهند بود. 


تاثیر ۱ ۳ اترای 3 کل واععه ی ۳ 


تا فانزم ۱,۱ نیز بصورت روابط ۳ فرض می‌شود: 
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مزا کر کی سروک فر زنل شلات هر طوقال اف رطف یج تباید ساخیان سین ۱۳۷ 


[ک501 1ک 501 زناس 10 
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۵( 
در اين روابط ۲ شعاع ذرم ٩و‏ و 650 کشش‌های 
سطحی آلیاژهای مذاب و جامد. ۷۳۹ و ۷*۶ حجم‌های مولی 
آلیازهای: مذاب.و امد تشد هه تروق ۵9 
کثش‌های سطحی جامدهای خالص ۸ و ظ و ۷ و 0و۷ 
نیز حجم‌های مولی آنها بوده که مقادیر آنها در مراجع موجود 

است. 
حجم‌های مولی آلیاژهای مذاب و جامد ۷14 و 750 


به صورت روابط زیر فرض می‌شوند: 
۱۰ و + رز 9نا7ا 
)۱۱( + ۷۵ -< ۷۶۵ 


کشش سطحی آلیاژهای مذاب و جامد. (0۶ وذل< 
01) . نیز طبق رابطه بوتلر (6 .1.۸.۷) بصورت زیر 
قابل محاسبه می‌باشد. 
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1 :0( لب ی ِ ِ 
۱۳( 
در این روابط 0 و یه کشش سطحی اجزای مذاب 
خالص, و 68۳7 و 16 کسر مولی عنصر 8 به ترتیب در 
سطح و توده می‌باشد که در این تحقیق برای محاسبه نمودار 
فازی یکسان فرض شدهاند. »۸ نیز سطح مولی جز خالص 
می‌باشد که طبق رابطه زیر از حجم مولی مذاب 0 ۷ 
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)۷/2 < بر 


۱( 
در این رابطه ول[ عدد آوو گادرو 6 یکی از اجزای ۸ یا 

ی ال اس شور سل تاد زو 
)۵7۳۳ نیز ارژی‌های گییس اضافی مولار 
جزئنی ۲ در سطح و توده می‌باشند که به صورت تابعی از 1 
و "108۳ یا 58 در نظر گرفته می‌شوند. انرژی گیبس اضافی 
مولار جزئی نوده مستقیماً با استفاده از روابط زیر به دست 


و عوهعهم - 31 | 
و وا( 6 1 0 
(۱۵( 
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و0 
)۱۹( 


انرژی گیبس اضافی مولار جزئی سطح نیز به صورت 
رابطه زیر به دست می‌آید: 


۷ 6 ار .۳1( 9 09 بت ا 09 


که برای آلیاژهای فلزی مذاب 0.83 < 8۳۴ و برای 
آلیاژهای جامد 0.75 < *8۳1 می‌باشد. به این ترتیب روابط 
(۱۲) و (۱۳). معادلاتی خواهند بود که با حل آنها می‌توان 
قفکار سرا تلبت وف 

در نهایت با جایگذاری پارامترهای فوق در معادلات 
(۸) و )٩(‏ می‌توان مقادیر انرژی آزاد سطح آلیاژهای جامد و 
مذاب را بدست آورد. با اضافه کردن این مقادیر به انرژی ازاد 
توده ماده (معادلات * و۷) انرژی آزاد کل محاسبه و در رسم 
نمودارهای نانو استفاده ی گرد [8]. 

در این مطالعه به منظور ساده کردن روابط فرضیات زیر 
قرمتحا بان تال یه اس 

کسر مولی عناصر تشکیل دهنده آلیاژ در سطح و توده 
یکسان فرض شده است. 

از وابستگی دمایی حجم مولی مذاب و جامد خالص 
صرف نظر شده است. 

کشش سطحی مذاب و جامد خالص مستقل از دما و به 


ترتیب برابر با مرتاه و بر؟1.250 در نظر گرفته شده است. 
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۱۳/۸ 


در ادامه دو مدل هایمینگ و تاناکا بر روی سیستم 
دوتایی ۸2-۳0 اعمال شده و نمودار فازی این سیستم برای 
نانوذره‌ای به قطر ۱۰ نانومتر و همچنین حالت توده آن رسم 
شبلة اشبرقا (شکل ۹ 


پسسسسته 9 200 


تساه ماتاهم‌نلا 


سس ۱۱۵۱۱۳۱۸۰۱۵۱ مهد( 
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شکل 1 نمودار فازی سیستم ۵۵-۳0 در حالت توده و نانو (با قطر ۱۰ 
نانومتر) بر اساس مدل هایمینگ و تاناکا 


با توجه به شکل (1) با کاهش اندازه ذره به سبب تاثیر 
فاکتور سطح, تقاط ذوب کاهش یافته و خطوط جامد و مذاب 
متسه مرک هگا سا ری این ها تا ی ی 
ناژ خر تاه هرفازی ربا تکرش سفق مخاول امه 
هم ری وه وی راهطا سا ع زو 
وی اب سکن ابفلا فوفیایی اش که یرای ساوه 
کرفن ووابط انم سفن آیه له ور ظر گ قسگاه اسگ 


نتیجه گیری 
در این تحقیق, الگوریتمی ارائه شد که با پیاده‌سازی آن در 


محیط نرم‌افزاری متلب. امکان رسم نمودارهای فازی برای 


تراچ 


رسم نمودارهای فازی ثانو برای سیستم‌های دوتایی در محیط نرم افزاری متلب 


سیستم‌های با قاعده و بی‌قاعده وجود دارد. اساس این 
لگوریتم پیدا کردن تابع انرژی آزاد مینیمم در هر دما بر 
حسب ترکیب شیمیایی و همچنین یافتن خطی مماس بر آن 
ست که در تمام بازه پایین‌تر از منحنی انرژی آزاد مینیمم 
قرار گرفته باشد. به اين ترتیب نقاط تقاطع این خط با منحنی 
نرژی آزاد مینیمم محاسبه و در ماتریسی ذخیره می‌شود. با 
رسم این ماتریس در تمام دماها و در طول بازه ترکیب 


شیمیایی می‌توان به نمودار فازی سیستم مورد نظر دست 


یافت. 

در ادامه با در نظر گرفتن معادلات متاثر از بخش نانو 
که از مها تبانگ :هاگ و تناها شایی شده و اضتال 
آنها در مقادیر ورودی نرم‌افزار نمودارهای فازی نانو مربوط 
به سیستم‌های دوتایی [06-5 و ۸۵۵-۳ رسم شد. نتایج نشان 
می‌دهد نمودازهای رسم شده بر اساس الگوریتم مورد استفاده 
در این تحقیق مطابقت خوبی با آنچه در مقالات این سه 
مسق اراک شیف ارت 

بعلاوه نتایج تایید می‌کند با دریافت مقادیر کوچکتر 
اندازه ذرات در ورودی این الگوریتم میزان افت خحط جامد 
و خط مذاب نیز بیش‌تر خواهد بود. بعبارتی خطوط جامد و 
مذاب در حالت نانو نسبت به حالت توده به سمت پایین 
حرکت می‌کند. زبرا پتابر اصول هندسی حاکم» کاهش اندازه 
ذره تا محدوده نانوی سبب افزايش نسبت سطح به حجم شده 
و تاثیرگذاری انرژی سطحی موجب کاهش نقطه ذوب می- 
گردد. علاوه بر اين با کاهش اندازه ذره در مقادیر ورودی نرم 
افزار در نمودارهای دوتایی ساده ملاحظه شد . محدوده 
دوفازی کوچک‌تر شده و مناطق محلول جامد نیز گسترش 
می‌یابد 
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